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Abstract. 

Although it has long been known that the proper quantum field theory description 
of critical percolation involves a logarithmic conformal field theory (LCFT) , no direct 
C ! consequence of this has been observed so far. Representing critical bond percolation 

■ as the Q — > 1 limit of the Q-state Potts model, and analyzing the underlying Sq 

. symmetry of the Potts spins, we identify a class of simple observables whose two-point 

functions scale logarithmically for Q — > 1. The logarithm originates from the mixing 
of the energy operator with a logarithmic partner that we identify as the field that 
creates two propagating clusters. In d = 2 dimensions this agrees with general LCFT 
results, and in particular the universal prefactor of the logarithm can be computed 



. exactly. We confirm its numerical value by extensive Monte-Carlo simulations. 

PACS numbers: 64.60.De 64.60.F- 05.50+q 11.25.Hf 
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1. Introduction 

The analysis of two-dimensional critical geometrical problems such as percolation or the 
self-avoiding walks (SAW) in terms of conformal field theory (CFT) involves features 
more complicated than those appearing in minimal (e.g., Ising or 3-state Potts) models. 
One might think that these complications are due to the non-local nature of the 
interesting geometric observables, such as connectivities of clusters and loops. This 
viewpoint is however somewhat misleading, since genuine non-locality would break the 
invariance under conformal transformations. In fact the above problems definitely are 
conformally invariant and their non- locality is only apparent: It can be traded for non- 
unitarity, by reformulating these models in terms of vertex models or supersymmetric 
spin chains. The origin of the difficulty in building a CFT description of percolation or 
SAWs is precisely this non-unit arity, which turns out to have deeper consequences than 
that of certain minimal models, such as the Lee- Yang edge singularity. 

One of these consequences is that the operator algebra underlying the lattice 
description — or the Virasoro algebra that emerges in the continuum limit — of percolation 
or SAW problems involves representations that are not fully reducible. This leads 
to indecomposability and, ultimately, to the appearance of logarithms in correlation 
functions. This feature follows from the fact that the models cited above have central 
charge c = 0, and hence it is necessary to go 'outside' their minimal Kac table (made of 
just the identity operator) to describe non-trivial physical observables. This extension 
leads to logarithmic CFTs (LCFTs) whose study has attracted considerable interest 
over the last few years, as their potential role in condensed matter as well as string 
theory applications has become more evident. 

The best known property of LCFTs at c = is the existence of a logarithmic partner 
to the stress-energy tensor, whose presence is necessary (under some circumstances at 
least) to avoid the 'c = catastrophe' pp. Associated with this partner is a universal 
number, the so-called '6 number' or logarithmic coupling. It was suggested that b is a sort 
of effective central charge [2] that could potentially be used to distinguish different c = 
theories [3]. While there has been a lot of work on the determination of the b number 
and its relation with abstract module properties of the Virasoro algebra [I], it is only 
very recently that methods have been devised to measure it in a numerical experiment. 
These methods [51 16] are however rather indirect — in particular, they require a proper 
quantization scheme — and thus to this day no simple physical procedure to determine 
b in a real experiment has been proposed. 

The indecomposability that leads to the existence of a logarithmic partner to the 
stress-energy tensor can be traced back to the peculiarities of the Q — > 1 (resp. n — > 0) 
limit that relate percolation (resp. SAW) to the Q-state Potts (resp. 0{n) spin) model. 
It is then possible to predict other striking consequences [TJ El E] that should be observed 
at c = 0, such as logarithmic terms in certain correlation functions. These logarithms 
arise from degeneracies, where several operators that are distinguishable at Q (resp. n) 
generic 'mix' in the limit Q — > 1 (resp. n — > 0). This mixing is the physical phenomenon 
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resulting from the indecomposability in the formal algebraic description. Nevertheless, a 
clear physical interpretation of such logarithmic correlators is often difficult to find, and 
also hard to access numerically since the logarithm multiplies a much stronger power 
law singularity. Direct consequences of the logarithmic nature of the underlying CFT 
have thus, to this day, remained unobserved. 

The purpose of this Letter is to present a new example of a physical observable in the 
Q-state Potts model which, one the one hand, exhibits a pure logarithmic dependency 
(with no multiplying power law) in the Q — > 1 percolation limit, whose universal 
prefactor we compute analytically in d — 2 dimensions, and, on the other hand, has 
a sufficiently simple geometric formulation in terms of percolation clusters to make 
possible a numerical study. And indeed our thorough Monte Carlo simulations nicely 
confirm both the logarithmic scaling and our result for the universal prefactor. 

2. Percolation, Potts model, and continuum limit 

We first recall the well-known reformulation of bond percolation as the Q — > 1 limit of 
the Potts model. The partition function of the Potts model reads 

z = e n ex p i K5 w) > c 1 ) 

where K is the coupling between spins o~i = 1,2,...,Q along the edges E of some 
lattice G. The Kronecker symbol S aij(Tj equals 1 if <jj = <7j, and otherwise. Universal 
properties depend only on the dimension d, and not on the precise choice of G. Although 
the main conclusions of this Letter should be valid in any dimension, for the sake of 
simplicity, we shall restrict ourselves to d = 2 in the following. The generalization of our 
results in higher dimensions will be discussed in the conclusion. We therefore suppose 
that G is the square lattice. The transition line — which gives rise to a critical theory 
for < Q < 4 — is then given by the selfduality criterion e K = 1 + y/Q. 

We can expand Z by rewriting the local Boltzmann weight as exp [K8 aua ^\ = 
1 + Yz^o-^Oj, with p — 1 — e~ K . The set of edges (ij) G E for which we have the 
term yz - ^,^ are called 'occupied bonds'. They define a graph H, whose connected 
components are known as Fortuin-Kasteleyn (FK) clusters [10]. Since spins belonging 
to the same FK cluster are aligned, we can perform the sum over {cr^} in ([I]) to obtain 

Z OC Q t)clustcrs(H) p ttedgCs(H) (1 - p yE\-iedgcs(H) _ £ 2 ) 

H 

In this formulation, the number of colors Q can be thought of as a real parameter, and 
the limit Q — > 1 yields a sum over bond percolation configurations, with a probability 
p per occupied bond. This model is critical for p = p c = |, and although the partition 
function Z — 1 is trivial, the correlation functions of the model capture the salient 
geometrical properties of critical percolation clusters. 

The continuum limit of the critical Potts model is described by a CFT. To establish 
this standard result, one first notices that the (outer and inner) hulls of the FK clusters 
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constitutes a gas of loops. This in turn defines a height model (the loops being contour 
lines of the height) which can be argued to renormalize towards a Coulomb Gas (CG), 
that is, a compactified free boson with Lagrangian density £ = ^r(V0) 2 , along with 
additional 'electric charges' at infinity. The stiffness g G [2, 4] is given by Q = 2+2 cos ^r, 
so that percolation has g — |. Over the past thirty years, this mapping has allowed 
physicists to compute many interesting geometrical properties for percolation (and for 
the Potts model in general), including crossing probabilities and critical exponents. More 
recently, the continuum limit has been studied by mathematicians under the name SLE K , 
where k — 16/ g. 

3. Symmetric group Sq and operators in the Potts model 

To understand how logarithms appear in percolation, the key idea is the study of the 
symmetric group Sq in the (formal) limit Q — > 1 (see also [7]). We expect scaling fields 
of the underlying CFT to transform as irreducible representations (irreps) under the Sq 
symmetry of the Potts model. 

Consider first observables acting on a single spin O(o~i), for Q integer. Obviously, 
any such operator can be decomposed onto a basis of Q generators as 0{ai) = 
Yla=i c a,3o-i,a- The action of Sq defines a representation of dimension Q, which however 
is reducible. Indeed, 5 aua can be decomposed onto an invariant 1 = Yla=x ^ it a, an d an 
irrep tp a = 8 aua — 1/Q with Yla=i fa, — of dimension Q — 1. Obviously 1 corresponds 
to the identity operator, while tp a is the magnetization, or order parameter operator, 
of the Potts model. The two-point function {tp a (o'i)<Pb(crj)) vanishes if % and j belong 
to different FK clusters (the sums over <7j and aj being independent). At the critical 
point, (Lp a (ai)(pb(o-j)) decays algebraically as r _2A ^, where the (bulk) critical exponent 
^■•fi = ( ' 6 ~ 9 gg g ~ 2 ^ can be computed within the CG setup [HI [12] . 

Nothing particular happens when one takes Q — > 1 in these expressions. However, 
the limit Q — > in (ipa(&i)<Pb(o'j)) is ill-defined, and although we will not discuss it in 
details here, this is actually responsible for the occurrence of logarithms at the level 
of the identity operator. The Q — > field theory of free (symplectic) fermions can 
be interpreted geometrically in terms of dense polymers, or spanning trees on G, and 
in this latter context several similar results were obtained using exact combinatorial 
methods [13l[TH[T5]. This theory is also related to the logarithmic form of the Gaussian 
propagator in d = 2 and to the asymptotic behavior of the equivalent resistance in an 
infinite network of resistors via the Kirchhoff theorem [16] . 

In the remainder of the Letter we focus on the Q — > 1 percolation case for which one 
needs to consider observables acting on two nearest-neighbors spins Oiji) = 0{p^ 0i+i) 
in order to recover logarithms. We impose the constraint o~{ ^ whence these 

observables are Q x Q (symmetric) matrices with zero elements on the diagonal. As 
before, the starting point is the basis elements S ai>a S ai+lt b which should be symmetrized 
in {o"j,o"j + i} in order to obtain symmetric matrices. The decomposition of this 
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representation of Sq is straightforward and we find 



w-^ .(i )9W - 1)9 m, (3 ) 



Symmetric matrices 0(<Tj ; <Tj-|_i) C(°'i; o 'i+l) = 
with zero diagonal 

where we denoted the representations by their dimensions. The explicit expression for 
the generators of these representations reads 

E(ai,a i+1 ) =S a ^ ai+1 , (4a) 
<p a (a h a i+ i) = S ai ^ ai+1 (<p a (vi) + <fa{cri+i)) , (4b) 

1pab(0-i, «7<+l) = S a ^ b (S ai> j ai+1>b + 5 CTi , 6 ^ +l!a - ^ a + ^ ~ Q(Q - \) ^) ' ^ 

where we recall that <p a {o~i) = b~ ai ,a ~ 1/Q- The 'scalar' E is obviously a one-dimensional 
irrep. The operator <p a satisfies ^2, a <p a = 0, and transforms like the 'vector' ip a . 
The 'tensor' ip a b = ipf, a satisfies V'ab — for any b. This last irrep thus has 

Q(Q — l)/2 — Q — Q(Q — 3)/2 independent elements, as expected from ([3]). 

The fact that we obtained three independent operators acting on two spins can 
be understood physically in terms of fusion of one-spin operators (see also [IT]). As 
(Tj 7^ <Tj, all three operators act on two distinct FK clusters. The operator with highest 
symmetry, ifj, will then impose that the two clusters propagate until they encounter 
another i/j operator. In section [5] we further develop this geometric interpretation. 



4. Percolation limit and Jordan cell for the energy operator 

We can already see from (j4j) that ip a b becomes ill-defined in the (formal) Q — > 1 limit. 
We shall now see that a well-defined limit is obtained by mixing ip a b with E. Let us first 
study these operators from a quantum field theory point of view. The energy operator 
is given by e(ri) = E(ai,ai + \) — (E), where we subtracted the bulk expectation value 
of E so as to obtain a well-defined scaling field. This subtraction does not change the 
representation theoretic considerations, so e(ri) remains an irrep. This field corresponds 
to the thermal perturbation of the Potts model, and its bulk scaling dimensions is 
A s = - — 1 [HI [12]. In the CFT with generic real Q, we expect the following form for 
the two point function 

(e(rM0)) = i(Q)(g-l)r- 2A ^), (5) 

where A(Q) is a regular function of Q, with a finite non-zero limit A(l) for Q — > 1. The 
reasons why (s(r)s(O)) should vanish at Q = 1 is very natural from a lattice point of 
view, and was recently argued in the context of bulk logarithmic CFT [TH] . 

The above fusion considerations imply that the 'tensor' operator ip a b can be 
identified, with the so-called 4-leg watermelon operator [12] (the two propagating 
clusters correspond to four propagating hulls). Its bulk scaling dimension A? = 
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follows from a CG computation [2D]. We thus deduce the form of the two- 



point function 

2A(Q) ( 1 

+ - ^ x r ~ 2A ^ Q \ (6) 

+ (Q-1)(Q-2)J ' W 

where A(Q) is again a regular function of Q when Q — >■ 1, and the factor 2/Q 2 is purely 



conventional. However, the Kronecker symbols combination is completely fixed by (jig) . 

In the formal limit Q — >• 1 the two-point function ([6]) diverges. To cure this, we 
introduce a new field 

4>ab(r) = Mr) + Q^rif a ^ b. (7) 

Its two-point function is easily computed and in order to have a finite Q — > 1 limit, we 
must require ^4(1) = ^4(1), and that A e = A^ at Q — 1. The latter condition is indeed 
satisfied, since for g = | the CG results read A £ = A^ = |. Assuming also the former, 
the two-point function of ip a b has a well-defined Q — > 1 limit: 



(<Mr)</U0)> = 2A(l)r~ 5 / 2 



4\/3 

(^ac + <Li + <^c + hd + $aJbd + <LAc) H log r 

7T 



(8) 



where we have used 



Q^l Q - 1 7T 

We have thus identified a logarithmic two-point function (j8| as a result of the mixing 
(JTj) of the energy operator e and the 4-leg operator ip of the Potts model. This mixing is 
consistent with recent algebraic results in bulk LCFTs [T9| |2T] , but was obtained here 
from very simple physical arguments based on the u Sq=i symmetry" of the theory. 

Logarithms in LCFTs can be associated with the non-diagonalizability of the scale 
transformation generator — i.e., the Hamiltonian in the usual radial quantization — of the 
theory. To show that this logarithm corresponds to a Jordan cell, we now study the 
change of the field ipaj, under a scale transformation r — > Ar. Recall [12] that e and ip^ 
are both primary operators for generic Q and thus transform as 

e(Ar) =A~ A£ £(r), (10a) 

Vu(Ar) = A- A ^ ab (r) . (10b) 

Using this, one can readily show that at Q = 1 

?MAr) = A" 5 / 4 ( <Mr) + ^logA e(r) ] . (11) 



The field ip a f, is therefore mixed with the energy operator e(r) after a scale 
transformation. We stress that the appearance of logarithms in the two-point function 
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(a) (b) (c) 

Figure 1. Percolation configurations contributing to (a) Po(»~)> (b) Pi(r) (one cluster 
propagating a distance r = \r\ — r%\), (c) P2<V) (two propagating clusters). 

OH]) of the field ip is fully compatible with scale and conformal invariance of the critical 
theory. The point is that some operators are not pure scaling operators, and get mixed 
with others under a scale transformation. In other words, the scale transformation 
generator (or Hamiltonian) is non-diagonalizable, with a rank-2 Jordan cell mixing the 
two fields ipab and e. 

5. Geometrical interpretation 

We now show that flHJ) has a very nice physical meaning in terms of geometrical 
observables in the percolation problem. Let i\ and + 1 be two nearest-neighbor points, 
separated from two other nearest-neighbor points 22 and 22 + 1 by a distance r = \r\ — r%\- 
To be consistent with the bulk calculation of section HJ we consider percolation defined 
on a torus {i.e., with doubly periodic boundary conditions) and we assume that r is much 
smaller than either period. We define P^ as the (position independent) probability that 
two neighboring points belong to different FK clusters. Let Po(^) be the probability that 
+ ii-, ^2+1} belong to four different FK clusters. Note that Po(r) is an increasing 
function of r, with limit (P^) 2 for r — > 00. Further let Pi(r) be the probability that 
{21,21 + 1,22,22 + 1} belong to three different FK clusters, of which one 'propagating' 
cluster contains one point from {ii, %\ + 1} and one point from {22, 22 + 1} — there are four 
ways to make this choice. Finally, let P2<V) be the probability that {ii, i\ + 1, 22, 22 + 1} 
belong to two different FK clusters, i.e., with two distinct clusters each containing one 
point from {ix,i\ + 1} and one point from {z 2 , «2 + 1} — this can be done in two ways. 
Configurations contributing to these probabilities are shown schematically in Fig. [TJ 

Using these probabilities, one can interpret geometrically the Potts operators 
previously introduced. For example, it is straightforward to show that (E) = ^pP^. 
Recall that the energy operator e was constructed by subtracting this expectation value 
from the operator E. A careful analysis shows that the two-point correlator of ip a b is 
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proportional to F 2 (r) 

2 / 1 

(ipabio-n, 0-i 1+ l)lpcd{0-i 2 , 0- i2+1 )) = — \ S ac Sbd + S ad 8 bc - - - (8 ac + 5 ad + 5 bc + 8 bd ) 



+ {Q - 1){Q - 2) ) >< ™ ■ < 12 > 

We infer from ([2]) that D^C?") ~ A(Q)r _2A '^ . Other correlation functions follow from 
the lattice description in the same way. A direct computation shows that (ip a b) = (0a) = 
0. This also follows from the representation theory of Sq, which moreover implies 
the vanishing of the 'crossed' correlation functions: (eipab) = (z(j>a) — (^ab'Pc) = 0. 
The vanishing of one-point functions and crossed two-point functions is consistent with 
fundamental CFT results. 

To analyze the Q — > 1 limit from a CFT perspective, we studied the correlation 
functions of the field ip a b defined in (J2J). One can repeat the very same steps from the 
lattice perspective. It is convenient to write ip a b( r i) = ^W^ij °"i+i) + q(q~ = 
' l l } ab{ (J i, 'i+i) — where we have introduced 

l/>db(<Ti, Ci+l) = ^,0^+1,6 + 8 a . jb 5 a . +ua - — - - (0 o (cr;, a i+ i) + 6 (cTj, cr i+ i)) . (13) 

Note that ip ab is not a scaling field (since (if>) = ^P^ ^ 0), whereas ip ab is. One can 
show that the two-point function of ip ab reads 

{iMn)<M>-2)} = 4 (Po(r) + Pi(r)) + P 2 (r) 



2 

(<$ac + S ad + he + Sbd) + 2(8 ac 5 bd + 5a Ac) 



_Q(Q-2) 2-Q 
whereas that of the corresponding scaling field i/j ab is, in the limit Q = 1 



, (14) 



(4>ab{ri)$cd{r2)) = 2 (5 ac + 5 ad + S bc + 5 M + 5 ac 5 bd + 5 a Ac) P 2 (r) 

+ 4[P (r-)+Pi(r)-2P 2 (r)-Py . (15) 

Comparing with (jSJ) we deduce that P2(r) ~ A(l)r~ 5 ^ 2 , as was of course expected 
from its relation to the 4-leg operator. Meanwhile, the logarithmic term in can be 
identified with 

P (r) + Px(r) - Pj ~ (e + ^ logr j r" 5/2 , (16) 

where we have added a subdominant non-universal (i.e., lattice dependent) term 9. 
Finally, the following combination 

P ( r ) + P l(r )-p2 2V3 
W = FW ~ # + -- logr, (17) 
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log r 

Figure 2. Monte Carlo estimation of the function F(r) defined in cq. (|17|) . Results 
are shown for 200 x 200 and 300 x 300 square lattices, with no significant difference. 

cancels out the dominant power law (r~ 5 / 2 ), leaving a pure logarithmic scaling which 
should be observable in numerical simulations (see below). The number ~ 1.1026 
is a universal constant that can be traced back to ((9]). Although the combination ( jTTl) 
may look slightly complicated, it is important to keep in mind that the logarithmic term 
we are after resides in the disconnected peat Po(r) — a similar observation holds true for 
LCFTs with other values of Q, such as Q — > 0. 

We have checked the validity of (TTTj) by performing extensive Monte Carlo 
simulations on square lattices of various sizes ranging from 150 x 150 to 300 x 300, with 
doubly periodic boundary conditions. We checked that different pseudo-random number 
generators — including a Mersenne-Twister algorithm [22] — led to consistent results. 
Statistics were obtained on ~ 10 3 independent runs of 10 7 percolation configurations 
each. Results are shown in Fig. [21 and are in good agreement with ( TTTj) . Careful 
extrapolations removing successively the first few short-distance points yields a slope 
1.15 ± 0.05 in good agreement with our prediction ~ 1.1026. 

Note that although all the calculations of this Letter were made in the bulk, the 
derivation in the boundary case presents only minor differences. In this case the scaling 
dimension of the energy operator should be replaced by A £ = 2, since e becomes 
degenerate with the stress-energy tensor T [23]. Similarly, the 4— leg watermelon 
exponent should be changed [24] to = — 2, and the operator ip ab would be 
proportional to the well-known logarithmic partner t(z) of the stress-energy tensor 
introduced by Gurarie [2] and Gurarie and Ludwig [3], [25] in their work on CFTs with 
central charge c = 0. However, it is easy to see that the limit ([9]) remains unchanged so 
we expect (fT7|) to hold true also if the points lie at a boundary. It would be interesting 
to check this numerically as well. 
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6. Conclusion 

We have found a simple geometrical observable in percolation that provides a lattice 
version of a logarithmic two-point function in a c = (L)CFT. The method we 
used is very similar to what was done by Cardy for disordered systems and for the 
0(n) model [7], and is far more general in that respect than the specific logarithmic 
solutions [131 EJ EES] found in the case of the free-fermion dense polymers case (Q = 0). 
Meanwhile, our result ( I17j) for percolation allows for a direct numerical verification, and 
the prefactor in front of the logarithmic term turns out to be universal. It should be 
noticed that this coefficient is actually closely related to what is known as logarithmic 
coupling or indecomposability parameter in the context of LCFT [31 [251 HI E] • We also 
note that, remarkably, in all the other examples we have studied, logarithmic terms 
tend to appear in disconnected observables such as Fo(r). This is similar to results for 
disordered systems [TJ, and we will get back to this issue elsewhere. 

The main parts of our derivation — notably the representation theory of section [3]— 
are not restricted to the d = 2 dimensional case. Indeed, we expect A? = A e at Q = 1 
also for d > 20, and only the derivative (J9]) will change. Accordingly we expect (JTTt) to 
remain correct for d > 2, albeit with a different universal prefactor that could be, in 
principle, computed in a e = 6 — d expansion. 

Our results are not restricted to percolation. In particular, studying the Q — ^2 
limit should also yield logarithmic observables, involving geometrical properties of Ising 
spin clusters. Higher-rank correlation functions and observables acting on more spins — 
including the generalization of 2n-leg watermelon operators to d > 2 — can be worked 
out along the same lines. There are indications that matching these more physical 
observations with formal algebraic developments should lead to further progress in our 
understanding of LCFTs. 

In conclusion, it is important to stress that logarithmic terms such as those we 
have identified would not be present for generic Q, and occur solely because of the 
special degeneracies present at Q — 1. This is of course quite different from logarithmic 
dependencies in other non-local quantities — see e.g. [27J [28] — which are obtained as 
derivatives of correlation functions with respect to the Boltzmann weights (such as Q). 
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% This implies that the fractal dimension g?rb of the so-called "red bonds" (also called "cutting bonds" ) 
is related to the thermal exponent v via ^rb = v . This is indeed a well-known percolation result [26] , 
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